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This paper presents the detailed modeling of various components of a grid 
connected hybrid energy system (HES) consisting of a photovoltaic (PV) 
system, a solid oxide fuel cell (SOFC), an electrolyzer and a hydrogen 
storage tank with a power flow controller. Also, a valve controlled by the 
proposed controller decides how much amount of fuel is consumed by fuel 
cell according to the load demand. In this paper fuel cell is used instead of 
battery bank because fuel cell is free from pollution. The control and power 
management strategies are also developed. When the PV power is sufficient 
then it can fulfill the load demand as well as feeds the extra power to the 
electrolyzer. By using the electrolyzer, the hydrogen is generated from the 
water and stored in storage tank and this hydrogen act as a fuel to SOFC. If 
the availability of the power from the PV system cannot fulfill the load 
demand, then the fuel cell fulfills the required load demand. The SOFC takes 


required amount of hydrogen as fuel, which is controlled by the PID 
controller through a valve. Effectiveness of this technology is verified by the 
help of computer simulations in MATLAB/SIMULINK environment under 
various loading conditions and promising results are obtained. 
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1. INTRODUCTION 

Due to the considerations of environment condition, applications of renewable energy sources 
(RESs) are more promising; that makes environment pollution free. Also, their availability is cost free and 
continuous [1], [2]. Numerous RESs consist of photovoltaic system (PV), wind turbine system (WT) and 
micro-turbines, etc. which are considered as components of hybrid energy systems in the literature and also 
demonstrate applications of micro-grid [3], [4]. Due to the various seasonal and bad weather conditions such 
as temperature, wind speed, solar radiation and also geographical conditions, these structures are not worked 
properly. So, the solutions must be needed and find out. Hence, energy storage systems (ESSs) are suitable 
for the solution of the mitigation of wind effects, solar radiation fluctuations and also ESSs uphold the power 
and energy balance. Power quality also improves due to ESSs. Due to the fast variations of power, ESSs must 
contain a high power density as well as high energy density. So, it is required to keep more than one storage 
system for a hybrid energy storage system (HESS) [5]-[7]. The ESSs and battery banks (BBs) are efficiently 
used in hybrid energy systems (HESs). But, lifetime of batteries decreases due to the charging and 
discharging cycles [8], [9]. A secondary energy sources is required to enhance the supply energy reliability of 
HESs. Hence, a fuel cell (FC) is required to combine with the electrolyzer by giving a continuous supply to 
the load [10], [11]. The strategies of energy management consist of combination of PV, WT and FC 
comprising with electrolyzer as well as battery storage. These are most effective practice for quality of higher 
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energy requirement. Further, micro-grid comes into interface with the distributed generation and for the 
decentralized management, the grid is connected with this micro-grid [12]. Now a day’s DC loads are 
increasing in electric vehicles, residential, commercial and industrial buildings. So, the loads in power system 
might be DC dominated in the subsequent time. If AC grid supplies these type of DC loads, then AC/DC as 
well as DC/DC converters are required in near future. These converter equipments combine with the 
individual DC or AC grid, so the cost is increased as well as cause the additional losses. In the meantime the 
overall efficiency is reduced for the particular system. Numerous studies and configurations have been 
suggested in [13]-[16] for the application of the hybrid energy system. A HES normally consists of individual 
DC grid or AC grid that can beat the fact that DC loads connect with the DC grid [17] and AC loads connect 
with the AC grid, also DC or AC grids are associated with the bidirectional converter [18]. 

From [18], a hybrid system includes PV system, WT system as well as FC. So, the hybrid system 
reduced the power electronics equipments such as AC/DC or DC/AC converters in the separate DC or AC 
grid applications. Also, FC is operated separately without batteries. But, robustness, vector control, 
advantages of WT-DFIG are not suggested, also their operations modes are not considered and further pitch 
angle control is not considered. Amongst other, Tamalouzt et al. [19] have discussed a direct torque and 
reactive power control of grid connected doubly fed induction generator for the wind energy conversion. 
Tazerart et al. [20] have described a direct torque control implementation method with losses minimization of 
induction motor for electric vehicle applications with high operating life of the battery. Bhende et al. [21] 
have designed a permanent magnet synchronous generator-based standalone wind energy supply system. 
Onar et al. [22] have described dynamic modeling, design and simulation of a wind/fuel cell/ultra-capacitor- 
based hybrid power generation system. In HES though wind integration looks attractive in terms of 
requirement of active power still the cost and control technique are the major drawbacks of this technology in 
addition to noisy environment. 

The major contribution of this paper is modeling of a new configuration of grid connected HES 
consisting of PV, electrolyzer, storage tank and SOFC which are simulated together. Particularly, the detailed 
modeling of hydrogen flow from electrolyzer to SOFC through storage tank which is not available in the 
literature has also been provided in this paper. Advantage of this technology shows the accessibility of grid 
(AC and DC). Hence, hybrid system increases power management capability of the grid. This paper provides 
the overall description of the proposed methodology. Further, detailed modeling and control strategy of each 
equipment of the system is described. A supervisory control strategy and also a management strategy of the 
system are presented. Finally, the effectiveness of the proposed HES is confirmed through the simulation 
under several situations such as variation in power demand and various meteorological conditions. 


2. MODELING OF COMPONENTS OF PROPOSED GRID CONNECTED HES 

The proposed model is shown in Figure 1. In this model the PV system is regulated by the 
maximum power point tracking algorithm (MPPT) and it is also interfaced with the boost converter to boost 
up the DC voltage. Then this DC voltage (Vpv_boost) is applied to DC/AC converter to convert DC voltage 
to AC voltage. Solar energy depends upon the weather conditions; so when the solar radiation is weak or it is 
absence for sometimes, then the SOFC comes into interface because integration of the system with SOFC 
makes it more sustainable. Further, when power generation is more from Sun for a long period, then the 
electrolyzer takes part to consume the excess electric power and generate hydrogen. This hydrogen is stored 
in storage tank in compressed manner and further utilized as a fuel of the SOFC. Therefore, the power 
management strategy of this HES is developed to minimize the effect of daily as well as seasonal variations 
due to the conditions containing climatic and geographical. So, RES results a polished output power, also 
satisfy the power requirement. Operating performances of HES depend on stability of the voltage of DC bus 
when disturbances/faults occur or when power demand fluctuation is rapid. When, voltage fluctuations is 
more, then power converters show more power losses and more injection of harmonics towards the grid. 
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Figure 1. Hybrid system connected to grid (Proposed Model) 


2.1. Modeling of photovoltaic (PV) system 

Mathematical model of the PV system is shown in [15], [16]. The power characteristic of the PV 
system is represented in [2]. To draw out the maximum power which is available in PV array, it is useful to 
operate the PV system at MPPT. Maximum power point tracking (MPPT) is a process where PV array, 
inverters which are connected to grid and other similar devices are employed to extract the maximum amount 
power. Hence, the PV module is designed with the help up MATLAB/SIMULINK by considering the 
following Equation (1). Also Table 1 contains specifications of the photovoltaic (PV) array. 
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Parameters considered to design the PV system such as Ip is the reverse saturation current of PV cell 
in [A], Isc is the short-circuit PV cell current in [A], Ipy/I,, is the output current of PV cell in [A], k is the 
boltzmann’s constant in [J/°K], a is the completion or ideality factor, q is the electron charge in [C], Rpis the 
PV cell containing parallel resistance in [Q], Rs is the PV cell containing series resistance in [Q], Ns is the 
no. of series cells in a string of the PV cell, Np is the No. of parallel strings, T is the temperature of the PV 
cell in [K], Vpy is the PV cell terminal voltage in volt in [V],Vyp is the voltage related to maximum power of 
the PV cell in [V], Vocis the PV cell open-circuit voltage in volt in [V]. 


Table | Specifications of Photovoltaic (PV) array 


Specification Value 
No. of series connected PV cell 96 

No. of Modules in parallel 66 

No. of Modules in series 05 
Open circuit voltage (Voc) 64.2 V 
Short circuit current (Isc) 5.96A 
Maximum Power of PV cell 305W 
Temperature (T) 25 
Irradiance 1000 


2.2. Modeling of electrolyzer 

Electrolyzer is used to decompose the water (H20) into two elements, first one is hydrogen and 
another one is oxygen by circulating the electric current in the electrolyzer containing two separate 
electrodes. The process is called electrolytic process or electrolysis [11]. 

The electrolysis equation of water is shown below: 


H,0(1) + electrical energy = H2(g) + ~ 02 (g) (2) 


Modeling, Control and Power Management Strategy of a Grid connected Hybrid ... (Sujit Kumar Bhuyan) 


1348 O ISSN: 2088-8708 


Parameters used for modeling of an electrolyzer are ie is the current of an electrolyzer in [A], ne is 
the no. of series connected electrolyzer cells, nH2is the amount of hydrogen generated in moles per second in 
[mol s '], F is the faraday constant in [C kmol "],nris the Faraday efficiency. 

Electrical current is directly related to the rate of production of hydrogen of the electrolyzer as 
shown in Faraday’s law given below [22]. 


— NFNcle 
nH, = Meee 3) 
The ratio of the actual amount of hydrogen generated by the electrolyzer to the theoretical value is 
represented as Faraday efficiency. 


0.09 755) 


ie B 


nr = 96.5el e a (4) 


From the Equations (3) and (4), an electrolyzer model is designed through the help of Simulink that is shown 
in Figure 2. 


Figure 2. A brief model of electrolyzer using Simulink 


2.3. Modeling of storage tank 

Hydrogen is produced by the electrolyzer, then this generated hydrogen acts as fuel of the SOFC. 
According to the output power, SOFC consumed the required amount of hydrogen. The difference between 
the generated hydrogen and required amount of hydrogen is stored in the storage tank in storage tank. 

In the hydrogen storage tank, the hydrogen is stored in compressed liquid or gaseous form. The 
hydrogen storage tank is designed according to the Equation (5) and tank pressure is directly calculated by 
considering the difference amount of hydrogen flow to the storage tank [22]. Figure 3 shows the brief model 
of hydrogen storage tank using Simulink. 


NH2RTp 


P, — Pi =Z 
b bi MH2Vp 


(5) 

Parameters used for modeling of the storage tank such as Z is the pressure compressibility factor, 
NH? is the amount of hydrogen stored in the storage tank in moles per second in [kmol s '], MHz is the 
hydrogen molar mass in [kg kmol '], Tb is the operating temperature in [°K], Ppi is the storage tank pressure 
at initial stage, Pp in storage tank pressure, R is the rydberg/Universal gas constant in [J (kmol °K) ~™'], Vp is 
the storage tank volume in [m°). 
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Figure 3. A brief model of hydrogen storage tank using Simulink 


2.4. Modeling of solid oxide fuel cell (SOFC) 

Parameters of the SOFC which are used to develop the mathematical model are given in [23], where 
B and C is the activation voltage constants are used to simulate in SOFC system in [V], CVis the conversion 
Hydrogen (kmoh) 

Methen (kmoh) 

is the Faraday’s constant in [C (kmol!’)], Ipc is the feedback current of SOFC in [A], k1 is the gain of PI, Kan 
is the anode valve constant in [./kmol kg(atm cS al | Kr is the modeling constant in [kmol SA”, Kuzo is 
molar constant of water valve in [kmol (atm s], Kuz is the molar constant of hydrogen valve in [kmol (atm 
s')], Koz is the molar constant of oxygen valve in [kmol (atm s')], Mrz is the hydrogen molar mass in [kg 
kmo!'], P20 is the water fractional pressure in [atm], Py2 is the hydrogen fractional pressure in [atm], Poz is 
the oxygen fractional pressure in [atm],Npo is the no. of fuel cells are in series in a stack, qu2 is the molar flow 
of hydrogen in [kmols”'], methane is the methane flow rate in [kmols”'], qoz is the input molar flow of oxygen 


in [kmols”], qi, is the input flow of hydrogen in [kmols”], qf: is the hydrogen output flow in [kmols”], 


qhz is the hydrogen flow amount that reacts in [kmols”], dis is the hydrogen flow amount that meets load 
change in [kmols”'], R is the Rydberg (Universal) gas constant in [(latm) (kmol K) 1, Rint internal resistance 
of SOFC in [Q], T is the absolute temperature in [K], U is the utilization rate,Van is the volume of the anode 
in [m], Veen is the SOFC dc output voltage in [V], T4 , T2 is the reformer time constants in [s],t3 is the PI 
controller time constant in [s], Ty is the hydrogen time constant in [s], Ty29 is the water time constant in [s], 
To is the oxygen time constant in [s], Zact is the activation over voltage in [V], Zonmic is the ohmic over 
voltage in [V]. 

The ratio between the hydrogen molar flow through the valve and partial pressure of hydrogen in the channel 
can be represented as [23]. 


factor in| |. Eo is the normal no load voltage in [V], E is the nernst immediate voltage in [V], F 


qH2 Kan 
4H2 — “an — K 6 
Py ¥Mn2 (6) 


Three important factors can be considered due to the flow of hydrogen: input flow and output flow of 


hydrogen and flow of hydrogen at the time of reaction. So, the equation can be developed by considering 
these factors. 


d RT i 
iP = p (dita — a2" — diz) (7) 


Reacted hydrogen flow rate is determined by considering the relationship between SOFC current and 
hydrogen flow. 


qiz = “SFE = 2K, Ipc (8) 
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The hydrogen partial pressure is determined by considering the laplace transform of the Equations (6) 
and (8). 


1 
/ : 
Pro = yet (atta — 2K lec) (9) 
where, 
Van 
Thi = nRT (10) 


In this way oxygen partial pressure as well as water partial pressure are determined. The summation of the 
activation over voltage, Nernst’s voltage, and the ohmic over voltage gives the result of SOFC polarization 
curve by considering the oxygen concentration and temperature to be constant. The voltage output of SOFC 
power plant can be indicated as: 


Veeu = E + Nact + Nonmic (11) 
where, 

Nact = —Bln(Clrc) (12) 
and 

Nonmic = —R™ Irc (13) 


So, Nernst instantaneous voltage can be obtained as: 


a RT Pu2vPo2 
E = No Ç + tog [aee] (14) 


The power plant of SOFC takes the hydrogen as a fuel to fulfill the requirement of power demand but the 
hydrogen is continuously supplied by the reformer for the SOFC stack operation. The model of the reformer 
can be represented as: 


qH2 = CV (15) 


methanol T1T2S?+(T1+T2)S+1 


A PID controller is used to control the flow of hydrogen according to the requirement of output power of 
SOFC system. SOFC output current is taken as an input of the SOFC to obtain the feedback control and is 
given by: 


req _ Nolrc 
m =e (16) 


The availability of the hydrogen of the reformer can be utilized to maintain methane flow by the help of PI 
controller. It can be expressed as: 


Amethane = (k + 2) (“alee za qi) (17) 


T3S 2FU 


A 5 KW SOFC is designed by using the above parameters as shown in Table 2 and Figure 4 shows a 
brief model of SOFC system using simulink [24]. In this proposed model three 5 KW SOFCs are used. 
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Table 2. Specifications of SOFC to design the mathematical model 


Specification Value 
Activation voltage constant (B) 0.04777 [A] 
Activation voltage constant (C) 0.0136 [V] 
Conversion factor (CV) 2 


Faraday’s constant (F) 


96 484 600 [Ckmol'] 


FC system internal resistance (Rin) 0.26664 [Q] 

FC absolute temperature (T) 343 [K] 

Hydrogen time constant (Ty>) 3.37 [s] 

Hydrogen valve constant (Km) 4.22 x10” [kmol (s atm)"] 
Hydrogen—oxygen flow ratio (rmo) 1.168 

Kr constant =No/4F 2.2802 x 10” [kmol (s A)! 
Line reactance (X) 0:0809 [Q] 

Methane reference signal (Qmethrer ) 0.000015 [kmol s$] 

No load voltage (Eo) 1.18 [V] 

Number of cells (No) 88 

Number of stacks (Ns) 1 

Oxygen time constant (Tg2) 6.74 [s] 

Oxygen valve constant, (ko2) 2.11 x10° [kmol (s atm)'] 
PI gain constant (k1) 0.25 

Reformer time constant (T4 , T2, T3) 15 [s] 

Universal gas constant (R) 8314.47 [J kmolK`'] 
Utilization factor (U) 0.8 

Water time constant (T420) 18.418 [s] 


Water valve constant (Kmo) 


7.716 x 10° [kmol (s atmy*] 
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Figure 4. A brief model of SOFC system using simulink 


3. CONTROL AND POWER MANAGEMENT STRATEGY 

Figure 5 shows the detailed block diagram representation of proposed HES as shown in Figure 1. 
Figure 5(a) shows the block diagram of HES, where maximum power can be tracked from PV system with 
the help of MPPT. This power is utilized to fulfill the load demand and the extra power of the PV system is 
given to the electrolyzer, which is utilized to extract the hydrogen. This hydrogen further used as a fuel of the 
SOFC. The PV system with MPPT is connected to DC-DC converter (boost converter) to boost up the 
voltage and the boost converter is triggered by the gate pulse generated by the MPPT. Again SOFC is also 
connected to another boost converter to boost up the voltage of SOFC as shown in Figure 5b. The purpose of 
voltage source converter (VSC) is used to synchronize the system voltage with the grid voltage. So, the grid 
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voltage and frequency are used for the operation of DC-AC converter (Inverter). VSC consists of two control 
cascade loops, which are used for regulation of DC bus voltage where, first one is the voltage controller of 
DC-link and last one is the internal current loop. 


3.1. Voltage controller of DC-link 

The DC-link voltage is used to control the active power with its desired value. Output of the DC- 
link PI controller (ig*) acts as a reference of the current PI controller (active current controller) that is shown 
in Figure S(c). 


DCUNKPL yt = + cuRRENT | i 
cowo 7 > CONTROLER Yy ie 
a 
enerator |” 
wy i | CURRENTPI | + i 
it > CONTROLLER ge f 
7 
sa 
(a) Block diagram of Hybrid Energy System (HES) (c) Block diagram representation of Voltage Source 
Inverter (VSC) 
-sofe boost 
Controlling 
Syal 
Required 
Amount fe —> 
given to SOFC 
(b) Block diagram representation of Hydrogen flow (d) Block diagram representation of Fuel Cell 
from electrolyzer to SOFC Controller 


Figure 5. A detailed block diagram of Hybrid Energy System (HES) 


3.2. Internal current loop 

The d-—g control (synchronous reference frame) is used for transformation of grid voltage (V, Vp, 
V.) and grid current (i+ i,, ic) to d-q reference frame which rotates synchronously along with grid voltage. 
With the help of d-q control, the DC values are obtained by controlling the variables. So, the DC values are 
easier for design and control. Grid voltage derives the phase angle by the help of phase locked loop (PLL) 
which can be used for abc — dq transformation for synchronization of grid voltage with controlled current. 
PLL is a technique which can be used to derive the phase angle from the grid voltage. From Figure 5(c), it 
can be seen that a zero value is set to the reference reactive current (i,*). So, power factor is maintained to be 
one (unity). The current PI controllers provide the voltage outputs (Vq*, V,*). These voltage outputs act as 
input of the PWM generator to produce the pulses which is used to trigger the DC-AC converter [25]. 

From the Figure 5(b), it can be seen that buck converter is connected across the electrolyzer. As 
described before the extra power of the PV system is given to the electrolyzer, where water is divided into 
hydrogen and oxygen and this hydrogen acts as fuel of SOFC. The requirement of the amount of hydrogen 
for the SOFC depends upon the load demand and the difference amount of hydrogen (hydrogen generated by 


Int J Elec & Comp Eng, Vol. 8, No. 3, June 2018 : 1345 — 1356 


Int J Elec & Comp Eng ISSN: 2088-8708 O 1353 


electrolyzer — hydrogen required by SOFC) is stored in storage tank. Figure 5(d) shows the block diagram of 
fuel cell controller, where the voltage signals of boost converter (Vsofc_boost) is given to the PID controller. 
Then, PID controller regulates the valve according to the load demand and decides how much amount of 
hydrogen flow to the SOFC through the valve. The values of PID controller are decided by self tuned mode 
(P=0.00013317, I=3.3925227e-05, D=-6.254254172870e-05). 


3.3. Operational control strategy of HES 

The neediness of the operational control strategy is to manage the relationship between the load 
demand and hybrid generation system. The operational control strategy follows four possible paths to 
manage this. 

a. When PV power (Ppv) is greater than load demand (Pload), then PV power is supplied to the load 
according to its requirement and the extra power of PV is supplied to the electrolyzer. By the help of this 
power electrolyzer divides the water into two halves (hydrogen, oxygen). Then, hydrogen is stored in the 
storage tank and later it is supplied to SOFC through the valve. 

b. When PV power (Ppv) is less than load demand (Pload), then PV power is not sufficient to fulfill the load 
demand. So, extra power is required to fulfill the load demand with the PV power and that extra power is 
supplied by the SOFC. 

c. When PV power (Ppv) is equal to the load demand (Pload), then total PV power is supplied to the load, 
no extra power of PV supplied to the electrolyzer. 

d. When PV power (Ppv) is equal to zero, then no PV power is supplied to the load and at that time SOFC 
fulfills the requirement of the load demand. 


4. SIMULATION RESULTS AND DISCUSSIONS 

The proposed hybrid energy system (PV/SOFC-Electrolyzer) has been designed under 
MATLAB/Simulink environment and also different simulations are carried out to verify the performances. 
Figure 6(a) shows the generation of hydrogen (H3) by the electrolyzer with the help of electrolysis process, 
where extra power of PV acts as the input of the electrolyzer after fulfilling the load demand and this 
hydrogen is stored in storage tank. Then, the hydrogen is utilized as a fuel of SOFC. Figure 6(b) shows the 
flow of hydrogen (H,) to the SOFC according to the load demand and it is controlled by the PID controller, 
which is shown in Figure 6(c). The obtained voltage from SOFC is boosted up by the help of DC/DC 
converter (boost converter) which is shown in Figure 6(d). The boosted voltage (Vsofc_boost) is used to 
fulfill the required demand, when the load demand is greater than PV power. 

Figure 7 presents the voltage and current wave forms of the hybrid energy system fed to the load. 
Figure 8 presents the voltage, current and power wave forms of 3-phase load having the capacity of active 
power (P) = 8000 W, inductive power (Q1) = 2000 VAR and capacitive power (Qc) =100 VAR. 

From the above simulation results, it is observed that H, flow to the SOFC and PID controller output 
track the load variation perfectly. In addition, the output voltage of boost converter connected to SOFC is 
matched to the required voltage for the load. Further, the output voltage and current of hybrid energy system 
and the voltage, current and power of 3-phase load are found to be satisfactory without having unwanted 
distortions. So, this proposed system can fulfill approximately different load requirements, while making the 
cost of this system cheaper which can be affordable by the consumers. This system is also environmental 
friendly and provides clean energy sources. 
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5. CONCLUSIONS 

This paper presents detailed modeling of a hybrid energy system consisting of PV/FC/Electrolyzer 
combined with a storage system. Different characteristics and system configuration of the various 
components are given. Also power management and the overall control strategy of the hybrid energy system 
are represented. The PV system is the main source of power generation, after fulfillment of the load demand 
by the PV generation the extra power of PV is utilized to produce H». The SOFC takes this H, as a fuel and 
acts as back up generation. Simulations are carried out under severe conditions and the system performances 
are verified as well as the results are obtained which show the performance of control strategy. In this 
proposed system, PID controller decides the amount of hydrogen flow to the SOFC through the valve 
according to the power demand. Also the hybrid energy system is connected to the grid with the help of VSC 
which results better performances of the whole system. Further, the results are verified and found to provide 
satisfactory performances. Hence, a better power management strategy is obtained to fulfill the load demand 
under the critical and unfavorable conditions. 
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